Aging Continuous Time Random Walks. 
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We investigate aging continuous time random walks (ACTRW), introduced by Monthus and 
Bouchaud [J. Phys. A 29, 3847 (1996)]. Statistical behaviors of the displacement of the ran- 
dom walker r = r(t) — r(0) in the time interval (0, t) are obtained, after aging the random walk in 
the time interval (— t a ,0). In ACTRW formalism, the Green function P(r,t a ,t) depends on the age 
of the random walk t a and the forward time t. We derive a generalized Montroll-Weiss equation, 
which yields an exact expression for the Fourier double Laplace transform of the ACTRW Green 
function. Asymptotic long times t a and t behaviors of the Green function are investigated in detail. 
In the limit of t 3> t a , we recover the standard non-equilibrium CTRW behaviors, while the impor- 
tant regimes t <C t a and t ~ t a exhibit interesting aging effects. Convergence of the ACTRW results 
towards the CTRW behavior, becomes extremely slow when the diffusion exponent becomes small. 
In the context of biased ACTRW, we investigate an aging Einstein relation. We briefly discuss aging 
in Scher-Montroll type of transport in disordered materials. 

PACS numbers: 05.45.-a, 74.40+k, 05.40.Fb 



INTRODUCTION 



Diffusion and relaxation in strongly disordered systems 
exhibits in many cases anomalous behaviors Jl], ||, [|. 
For example the diffusion of a test particle may become 
anomalous, namely the mean square displacement be- 
haves like (r 2 ) ~ t a and a ^ 1. A random walk frame- 
work, widely applied to describe anomalous diffusion is 
the continuous time random walk (CTRW) j§. CTRWs 
are used to model many physical and chemical processes, 
for example: charge transport in disordered systems 

transport in low dimen- 
anomalous diffusion in 
, the chemical reaction 
, and blinking behavior 




protein folding dynamics 
sional chaotic systems 
a metallic supercooled liquid 
of CO binding to myoglobin | 
of single quantum dots [^4[ |l5| 

Anomalous diffusion processes may exhibit aging |l6| , 
|l7| , |i~9| , |2f| , where vaguely speaking the age of the 
process controls the statistical properties of the random 
walk. Aging in diffusion processes yields an interesting 
insightful perspective on dynamics in disordered medium, 
and more generally is used as a tool to probe complex sys- 
tems such as spin glass, Anderson insulator, and colloidal 
suspensions (see |2jj for a brief review). Monthus and 



Bouchaud 1 22 introduced a CTRW framework, which ex- 
hibits aging behaviors. We call this generalized CTRW, 
aging continuous time random walk (ACTRW). In this 
paper, we investigate properties of biased and non-biased 
ACTRWs in detail, in particular an exact expression for 
the Fourier double Laplace transform of the Green func- 
tion is obtained, and asymptotic behaviors are investi- 



gated. 

ACTRW might yield a phenomenological description 
of aging self diffusion dynamics in glasses [E2| , below the 
glass transition temperature T c and then a = T/T c is 
temperature dependent (and see 23 |24|, |25|, [27) for 
related work). 



3 m 

showed that ACTRW de- 



Recently 

scribes dynamics of an intermittent chaotic system. More 
generally, we expect that any random walk described by 
a CTRW process will exhibit aging, provided that aging 
initial condition discussed in the manuscript are satis- 
fied. We note that other stochastic approaches to aging 
dynamics are based on a non-linear diffusion equation 
p9[ , and a generalized Langevin equation Q . 

This paper is organized as follows. In Sec. (|J) CTRW 
and ACTRW are introduced. We then derive an exact ex- 
pression for the Fourier double Laplace transform of the 
ACTRW Green function Pjr, t a ,t), thus generalizing the 
Montroll-Weiss equation p) to include the effect of the 
age of the random walk (see Sec. [ill]). Our result is based 
in part on recent work of Gordeche and Luck [ pT| who 
investigated renewal theory for processes with power law 
waiting time distributions. In Sec. IV we derive asymp- 
totic behaviors of the Green function P(r, t a , t) which are 
analyzed in detail. We then consider biased ACTRWs 
(Sec. ^) and discuss an aging Einstein relation. Finally 
we briefly consider aging in Scher-Montroll type of trans- 
port as a possible application of this work. Note that a 
small part of our results was reported in |^8| , in the con- 
text of aging in chaotic transport. 



II. CTRW AND ACTRW 
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One of the best well known random walk models is the 
CTRW introduced by Montroll and Weiss Q. It de- 
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scribes a large class of random walks, both normal and 
anomalous and can be described as follows. Suppose a 
particle performs a random walk in such a way that the 
individual jump r in space is governed by a probability 
density function (PDF) /(r), and that all jump vectors 
are independent and identically distributed. The char- 
acteristic function of the position of the particle relative 
to the origin after n jumps is /"(k), where /(k) is the 
Fourier transform of /(r). Unlike discrete time random 
walks, the CTRW describes a situation where the wait- 
ing time t between jumps is not a constant. Rather, 
the waiting time is governed by the PDF ip{t) and all 
waiting times are mutually independent and identically 
distributed. Thus, number of jumps n is a random vari- 
able. 

Let Pmw{ v , t) be the Green function of the CTRW, the 
Montroll- Weiss equation yields this function in Fourier- 
Laplace (k, u) space: 



Pmw (k, u) 



l-ip(u) 



1 



u l-/(k)^(«) 



(1) 



All along this work we will use the convention that the ar- 
guments in the parenthesis define the space we are work- 
ing in, thus tp(u) is the Laplace transform of ip(t). Prop- 
erties of Pmw (r, t) based on the Fourier-Laplace inver- 
sion of Eq. (|l|) are well investigated, see || |33| and Ref. 
therein. In particular, it is well known that asymptotic 
behavior of Pmw (r, t) depends on the long time behavior 
of tp(t). Two classes of processes are usually considered. 
The first is the case when all moments of i/j(t) are finite, 
the second class is the case where ijj(t) is moment-less, 
corresponding to a situation where oc and 
< a < 1. 

An important assumption made in the derivation of 
Eq. ([!]) is that the random walk begun at time t = 0. 
More precisely, it is assumed that the PDF of the first 
waiting time, i.e., the time elapsing between start of the 
process at t — and the first jump event is ip(t). Thus the 
Montroll- Weiss CTRW approach describes a particular 
choice of initial conditions, called non-equilibrium initial 
conditions. The limitation of CTRW theory to a very 
particular choice of initial conditions, was an issue for 
debates in the early seventies jiij . 

Monthus and Bouchaud Q introduced a CTRW for 
an ongoing process, where the random walk process is as- 
sumed to start at some time t = —t a , long before start of 
observation at time t = 0. In Fig. |l| a stochastic realiza- 
tion of number of jumps in such a process is shown. For 
such a random walk the Green function is denoted with 
P(r,t a ,t) and r is the displacement in the time interval 
(0, t). Using scaling analysis, (22j have investigated basic 
properties of this random walk, mainly the behavior of 
the Fourier transform of the Green function. For a < 1 
the random walk exhibits interesting aging effects, hence 
as mentioned we call it ACTRW. 




FIG 

-l 



1: Number of jumps i in a renewal process with ip(t) = 
i -1 ' 2 (l+t) _1 , i.e. a = 1/2. The process starts at t = —t a 



7r t ' (.i+tj , i.e. a - 
with t a = 4000. Observation of the process begins at time 
t — 0. ti is the time elapsing between t = and first jump 
event in the forward time interval (0, t), in the Fig. t\ — 2553. 



III. ACTRW:- GENERALIZED 
MONTROLL-WEISS EQUATION 

The ACTRW describes the following process, a parti- 
cle is trapped on the origin for time t%, it then jumps to 
ri, the particle is then trapped on ri for time t^, and 
then it jumps to a new location; the process is then re- 
newed. Thus, the ACTRW process is characterized by 
a set of waiting times {t\, t n , ...} and displacements 
{ri, • • • , r n • • •}. The time elapsing between start of ob- 
servation at t = 0, and the first jump event is denoted 
by t\. Here we denote the PDF of the first waiting time 
t\ with ht a (t\). In ACTRW the random walk started at 
t = —t a , before the start of observation at t = 0, there- 
fore h ta {t\) depends on age of the process t a . The waiting 
times with n > 1 are independent and identically 

distributed with a common probability density 4>(t). The 
jump length {ri, • • • , r n • • •}, are independent identically 
distributed random variables, described by the probabil- 
ity density /(r). 

In contrast, in the Montroll- Weiss non-equilibrium 
CTRW, the age of the process is zero t a = 0. And, for 
that case ht a {t%) = ^>(ti). 

Recently, Gordeche and Luck |3l| investigated statis- 
tical properties of fractal renewal processes, among other 
things they obtain h ta (ti). Let h s (u) be the double 
Laplace transform of ht a {t{) 



h s {u) 



dii / dt a h ta (t{) e 



-t a s—t\u 



(2) 
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then according to |3l| 



h s {u) 



1 - iP(s) 



(3) 



In Appendix A we re-derive Eq. (^|) using a method 
which slightly differs than the one used in |pl[ . 

Two types of behaviors are found for h ta (ti) The 
first case corresponds to a situation when average waiting 
time (t) = J °° tip{t)dt is finite, and then in the long aging 
time limit one obtains 



lim h ta (<i) 



g mat 

(t) ' 



(4) 



This type of initial condition is called equilibrium initial 
condition, it was investigated previously in the context 

Here we will 



of CTRW and related models 1 34 
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mainly consider the second case corresponding to a power 
law waiting time PDF 

cx t- {1+a) with < a < 1, (5) 

when t is long. In Laplace t — ► u space Eq. (||) reads 

i/j(u)~l-Au a , (6) 

where u is small, and A is a positive parameter |S7| . For 
example the one sided Levy PDFs whose Laplace pair 
is ip(u) = exp(— Au a ), or if>(u) = 1/(1 + Au a ) discussed 
below, belong to the class of functions described by Eqs. 
(||,|]). In the limit of long aging times, t a » A 1 / 01 , these 
kind of probability densities yield 



M*0 



sin(7ra) 



t? (tl+ta) 



(7) 



Note that this expression is independent of the exact form 
of ip{t), except for the exponent a. When a —> 1 the mass 
of the PDF h ta (t\) is concentrated in the vicinity of t\ — > 
0, as expected from a 'normal process'. Eq. (0) shows 
that as age of the process becomes older, we have to wait 
longer for first jumping event to occur. In a physical 
process, this may correspond to a particle in a disordered 
system which searches for a local energy minima in time 
interval (— t a , 0). In this case the longer the search takes 
place the deeper the minima found, hence in statistical 
sense the release time becomes longer as the process is 
older. In what follows, we will also use the double Laplace 
transform of Eq. (Q) : 



h s (u) 



s a (u — s) 



(8) 



This equation can be derived by inserting the small u 
and s expansion of tp(u) and ip(s) given in Eq. (^3j), in 
Eq. (1). 

Let P(r,t a ,t) be the Green function of the random 
walker, where as mentioned 



r = r(t)-r(Q) 



is the displacement in the time interval (0, t). Hence, 
clearly initially r = at time t = 0. Let (i) p n (t ai t) 
be the probability of making n steps in the time interval 
(0, t) and [ii) P(k, s, u) be the double-Laplace -Fourier 
transform (r — > k, t a — > s, t — » u) of P(r, t a , t). Then 



P{lL,S,u) = Y J Pn{s,u)f n (k) 
n=0 



(10) 



where p n {s,u) is the double Laplace transform of 
p n (t a ,t). As mentioned, / n (k) in Eq. ( p^| ) is the charac- 
teristic function of a random walk with exactly n steps. 
Using the convolution theorem of Laplace transform we 
obtain 



p„{s,u) = 



l — sh s (u) 



71 = 

n > 1. 



Hence inserting Eq. ( pd| ) in Eq. (|To|), using Eq. 
summing over n, we find the exact result 



(11) 



and 



P(k,s,u) 



1 

su 



[^{u)-^{s)\ [l-/(k)] 



tt (« - s) [1-iP (s)] [1 - ^ («) / (k)] ■ 

f(12) 
) is a generalization of the Montroll- Weiss equa- 
for ACTRW. Note that P (k = 0, s, u) = l/(su) 
as expected from the normalization condition. 

It is useful to rewrite Eq. (|l2|) in terms of the Montroll- 
Weiss Eq. (|l|), and Po(s, u) in the first line of Eq. (|ll|): 

P(k, s,u) =p (s,u) + h s (u)f(k)P M w(Ku). (13) 

The first term on the right hand side of this equation, de- 
scribes random walks where the particle does not leave 
the origin (i.e. n = 0). The second term describes 
random walks where number of steps is greater than 
zero, it is given in terms of a convolution of h s (u)f(k) 
with the Montroll- Weiss equation. This is expected 
since the only difference between ACTRW and the non- 
equilibrium CTRW, is the first waiting time distribution. 

If the process is Poissonian, tp(t) — exp(— t), the Green 
function P(r,t,t a ) is independent of the age of the ran- 
dom walk t a . To show this we insert 



ip(u) 



in Eq. (02) and find 



1 



1 



P(k,a,u) 



1 1 



(14) 



(15) 



su+l-/(k) 

Inverting to the double time domain 

P(k,t a ,t)^e^ 1 - f ^ t . (16) 

This result is independent of t a as expected from a 
Markovian process. Assume that /(k) = 1 — m^lkl^ 1 + • • ■ 
(9) for small values of k and /j < 2, implying that the 
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random walks is non biased. In the long time limit 
P(k,t a ,t) ~ exp(— m^\k\^t), and either a Levy behav- 
ior (/i < 2) or a Gaussian behavior (fi = 2) is found, as 
expected from the Gauss-Levy central limit theorem [B8| • 
In what follows, we investigate cases when this standard 
behavior does not hold. 



IV. ASYMPTOTIC BEHAVIORS 

Let us now consider basic properties of ACTRWs. 
While Eq. ([[2]) is valid for a large class of random walks, 
including Levy flights (/i < 2), we will assume that vari- 
ance of /(r) is finite (fi = 2). Special emphasis will be 
given to the case when ip{t) is moment-less a < 1, since 
this regime exhibits interesting aging behaviors. 



A. Mean Square Displacement 

By differentiating Eq. ([l2]) with respect to k and set- 
ting k = 0, we obtain the moments of the random walk 
in a standard way. Assuming a non biased symmetrical 
random walk, we obtain 



(r 2 M) = 



h s (u)m2 
u [1 — tp (it)] ' 



(17) 



where m,2 — / r 2 /(|r|)dr is assumed to be finite. We 
consider power law waiting time PDFs as in Eqs (|^||) , in 
the limit where both u and s are small, their ratio being 
arbitrary, we find 



(r 2 (*,«)) 



(u a — s a ) m2 
s a (u - s) Au 1+a ' 



(18) 



As shown below one can invert this equation exactly to 
the double time domain. However it is instructive to 
consider two limits first. If u <C s, corresponding to 
t > we have 



<r 2 (s,u)) 



u 1 a rri2 
As ' 



While for s <C u, corresponding to t a 3> t, we have 
Inverting Eq. (|l9|) and Eq. (|2C|), we obtain 



m 2 t 



(r 2 (t a ,t)) 



/H'(l+n) 



t > t a 



(19) 



(20) 



(21) 



This result is valid provided that both t, t a » A x l a . In 
the limit t 3> t a we recover standard behavior found in 
non-e quilibrium CTRW j§. In the aging regime, t <^it a 
we find an interesting behavior. Independent of the expo- 
nent a, the mean square displacement increases linearly 



with respect to the forward time t, as found in normal 
diffusion processes. In addition, the diffusion is slowed 
down as the age of the random walk t a is increased. This 
behavior is expected, due to statistically longer release 
times, from the initial position of the particle, as the age 
of the random walk is increased. 

We now consider a specific choice of waiting time dis- 
tribution: 

^(«) = — -a — » ( 22 ) 
1 + Au a 



corresponding to |T( 

m 



A 



.4 



(23) 



where E a , a (x) is the generalized Mittag-Lefner function 
p4j. Inserting Eq. (p2|) in Eq. (17) we have 



(s,u)) 



7712 u 



A s a (u~ s)u 



1+a ' 



(24) 



Hence for this choice of waiting times, Eq. (|8|) is exact 
and not limited to the asymptotic regime. Inverting to 
the time domain using Eq. (171) we find 



(*»,*)> 



7772 sin (ira) 



A 



r(i 



(25) 

where <g> is the Laplace convolution operator with respect 
to the forward time t. We rewrite Eq. (J25|) as 



(r 2 {t a ,t))=t 



777 2 



a ^ir(i 



sin (7ra) f t/ta (t/t a - y) a 
y a (l + y) 



dy. 

(26) 



The solution of the integral is readily obtained, we find 

7772 1 



(*<»»«)> 



A r(l + a) 



(27) 



The right hand side of Eq. ( p7| ) describes the long time 
t, long time t a , behavior of a large class of random walks 
with waiting time PDF satisfying ip(t) oc t~( 1+Q ) [i.e., 
since the right hand side of Eq. (|27]) is the double inverse 
Laplace transform of Eq. (Q)]. Note that if a = 1 in 
Eq. (^), the random walk does not exhibit aging. 



B. Green Function 

In this subsection we investigate asymptotic properties 
of the Green function P(r, t a , i), by considering the con- 
tinuum approximation of Eq. (|l2j) . This approximation 
is expected to work in the limit where both the forward 
time t and the aging time t a are long. A proof of the 

fVE], for the 



validity of this approach, is given in Sec. 
one dimensional ACTRW. We assume a symmetric ran- 
dom walk, hence for small |k|, the following expansion is 
valid: 



/(k) ~ 1 



2 |k| d ' 



(28) 
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where d is the dimensionality of the problem. We also use 
the small Laplace variable u expansion ip(u) ~ 1 — Au a 
in Eq. (Ell). Inserting these expansions in Eq. (|l2|), we 
obtain 



P(k,«,«) 



s Q u - su" (w a - s°) 
s^+^u - s) + s a {u - s) Au a + |k| 2 ^ 



2d 

For convenience, and without loss of generality, we choose 
now to work in units where A = 1 and 7712/ (2<i) = 1. 

Inverting Eq. ( p9|) to the double time (t a ,t) - real 
space r domain we find 



P(r,t a ,t) -JJ (t B ,*)«(r) + 



sin (jra) 



iWw (r, t) , (30) 



where ® in Eq. ( |30| ) is the Laplace convolution operator 
with respect to the forward time t, and in this limit 



Po(ta,t) 



sin (to) 



da; 



(31) 



In Eq. ( p0| ) PamwC 1 ")*) is the long time solution of the 
Montroll- Weiss equation, i.e., the Green function of the 
fractional diffusion equation 0, 



,a-l 



(32) 



where £ _1 .F _1 is the inverse Laplace u — ► f inverse 
Fourier k — > r operator. 

The Green function Eq. ( |30| ) , is a sum of two terms. 
The first term on the right hand side of Eq. (^) is a sin- 
gular term (i.e., the S(r) term). This term corresponds 
to random walks where number of steps in time interval 
(0, t) is zero. Unless t ^> t a , this term cannot be ne- 
glected, since without it the Green function in Eq. ( |30| ) is 
not normalized. Thus ACTRW exhibits a behavior differ- 
ent than ordinary CTRWs, where realizations of random 
walks where number of steps is zero do not contribute to 
the asymptotic behavior. 

In one dimension, we have 



PAM\v{r, t) 



t 



a 



/2 



t 



\2/a 



(33) 



where l a /2(t) is the one sided Levy stable PDF, whose 
Laplace pair is exp(— w Q / 2 ). Hence the Green function 
solution of the ACTRW is 



P(x,t a ,t)~ Po (t a ,t)S(x) + 



sin (tto) 



f °te)° 



'a/2 



t\x\ 



-(2/a) 



2 l/a 

(34) 



Scaling Eq. (§|) with r = t/t a and q = \x\/t a / 2 we 
obtain a scaled form of the Green function: 

P{x, t a ,t) ~ jj (t a , i) <^ (a;) + 



sin (to) /2 



1 



1+2/a 



^ r '7 A" Pi 1 A^ Q /2 

(t — t'J (1 + t — r'J 



2/a r J ' 



(35) 



Below we analyze this expression in some detail. 
In d space dimensions we have J42J 



Pa 



A1W 



(r, t) = a 7r 



tt2Q 
H 12 



Besides the d = 1 case, the Fox function solution P\\ 
is not generally tabulated, hence this solution is rather 
formal, though asymptotic behaviors of Eq. (|3^) are well 
investigated |33|, Q . A practical method of obtaining the 
solution of Pamw{t, t) Eq . ([56]), using the inverse Levy 
transform, is given in |33j. Using this method, we find 
the integral representation of the aging Green function 
in d dimension 

P(r,t a ,t)~ Po (t a ,t)S{r) + 

1 sin (to) t -ad/2 x 
t 2 d n d / 2+1 a 



ds 



U s l+l/a+d/2 a 



,1/q 



-q 2 /(4s) 



(1+r'Y 



(37) 

where q — r/t a ' 2 and r = t/t a . Similar to the one di- 
mensional case, this solution shows the precise relation 
between ACTRWs and Levy stable laws. 

To conclude, Eq. ( |3C| ) shows that the asymptotic so- 
lution of ACTRW is a sum of two terms. The first is 
a singular term, and the second is a convolution of the 
distribution of the first waiting time and the asymptotic 
Green function of the non-equilibrium CTRW. 



C. Graphic Examples 

In order to better understand the asymptotic behavior 
of the ACTRW, we perform numerical calculations to 
obtain the non-singular part of the Green function for 
different values of a in one dimension. The one sided 
Levy stable probability density in Eq.(^3|) was obtained 
using a numerical inverse Laplace transformation method 
El, El. The calculated PDF was then used to evaluate 



the convolution integral numerically to obtain the non- 
singular part of the Green function according to Eq.(35). 
For mathematical details on one sided Levy stable laws, 
see Appendix in |3^] , and references therein. 

In Figs. 0-pL we present the calculated non-singular 
part of the Green function for different a. Fig. || shows a 
three dimensional plot of the scaled non-singular Green 
function, P[q,r) ■ t a / 2 , versus scaling variables, r = t/t a 
and q = \x\/t a / 2 , for a — 1/2. A smooth transition from 
the aging behavior when t < 1 to the non-equilibrium 
CTRW behavior when r 3> 1 can be clearly seen in this 
figure, where the time ratio t = t/t a is changed con- 
tinuously. In addition, we observe a monotonic increase 
of the non-singular Green function as the scaled time r 
increases. Note that as the scaled r is increased the sin- 
gular term is decreasing, hence we may think of this aging 
process, as if the singular part of the Green function, is 
feeding the non-singular part. 

Figs. |^-|| show the scaled non-singular Green function 
versus q = \x\/t a / 2 for a = 1/6, a = 1/2, and a = 5/6, 
respectively. In each figure, the scaled Green functions 
at several different rs are shown. A few general features 
can be seen in these figures. First, the Green function 
is clearly non-Gaussian for all cases, as we expected. A 
comparison between the shape of the Green function for 
a = 1/6 in Fig. | and for a = 5/6 in Fig. || clearly 
demonstrates that the deviations from Gaussian behav- 
ior are stronger for smaller a. In the limit of a — > 1 
(not shown), we obtain a Gaussian Green function. Sec- 
ond, for t 3> t a (i.e. r » 1), we recover the usual non- 
equilibrium CTRW behavior shown as the solid curves in 
Figs. §-g Finally, as a becomes small, the convergence 
towards the non-aging behavior when t/t a — oo becomes 
extremely slow. For example, the case for a = 1/6 in Fig. 
|3| shows a significant deviation from the non-equilibrium 
CTRW behavior (solid line) when t/t a = 10000. As a 
result, the Green function exhibits aging behavior even 
when t/t a is large, and an ACTRW treatment for dy- 
namics in this regime will be essential. 



D. Proof of Asymptotic Behavior for One 
Dimension 




FIG. 2: Three dimensional plot of the scaled non-singular 
Green function versus r = t/t a and q — \x\/t a ^ 2 for a = 1/2 
and one dimension. Notice a smooth transition from aging 
behavior r << 1 to Montroll- Weiss CTRW behavior found 
when r >> 1. 




FIG. 3: The scaled non-singular Green function versus q = 
\x\/t a/2 for a = 1/6 and for five different r = t/t a : r = 0.1 
(dot-dash), r = 1 (dash), r = 10 (dot dash), r = 10000 
(dashed), and r = oo (solid). The solid curve is the asymp- 
totic behavior of the non-equilibrium CTRW. Notice the non- 
Gaussian shape of the Green function and the slow conver- 
gence towards the non-aging CTRW behavior (the solid line), 
compared with the cases a = 1/2 and a = 5/6 shown below. 



We now prove the validity of Eq. ( |30| ) using a method 
developed in |39) . The main idea is to show that moments 
of the ACTRW, are in the asymptotic limit described well 
by Eq. (|3^) . For simplicity we assume a one dimensional 
symmetric random walk. 

The moment generating function /(k) is expanded 

k 2 fc 4 k 6 

/(k) = 1 -m 2 y +TO 4 — - m 6 — + ■■ •. (38) 

Where rri2, etc are the moments of the jumps. In- 
serting this expansion in Eq. (|I^) we obtain the small k 
expansion of the ACTRW moment generating function: 

P(k, s,u) = --^-([l + Q(u)] mX 

su u y i 



- [1 + Sl(u)\ (m 4 + QVL{u)ml) — 

k 6 

[1 + fl(u)\ [m 6 ~ 30tt(u)m 2 m 4 + 90Q 2 (u)r,4] — + ■ ■ ■ 

(39) 

where 

fl(u) =ip(u)/[l-ijj(u)]. (40) 

The moments (x n (s,u)) of the ACTRW are defined in 
the usual way 

P(k,s,u) = jT(x n (s,u))^. (41) 

n=0 
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FIG. 4: Same as Fig. H for a = 1/2. 




are the irrelevant parameters in this problem. Inserting 
Eq. @ in Eq. @ we have 

"<M,.)~£^W(^f 2 . (46, 

Since we are interested in the limit where t and i Q are 
large, the ratio t/t a being arbitrary, the corresponding 
Laplace variables u and s must approach zero their ratio 
being arbitrary. Therefore h s (u) in Eq. ( fhsl) is given by 
its asymptotic form in Eq. (|^). Inserting this expression 
for h s (u) into Eq. (ff6|), setting m 2 = 1, and then sum- 
ming over n, we find an ex pre ssion for P(fc, u, s) that is 
the same as Eq. (29). Eq. (29) when transformed yields 
Eq. (^) . To conclude we showed that Eq. (|2^) describes 
the small s, u behavior of the ACTRW moments, hence 
it follows that Eq. (^) describes the long time t and 
t a behavior of the ACTRW Green function, the set of 
moments m^rriQ etc are unimportant in this limit. 



E. Behavior on the Origin 

Using Eq. ( |33|) we investigate the nonsingular part of 
the ACTRW on the origin. For d = 1, we find 



P(x,t a ,t)\ x=0 =t- a / 2 g(-), 



(47) 



FIG. 5: Same as Fig. H for a = 5/6. 



Comparing Eq. (B9h with Eq. (ttlh we have 



1 



<*°> = 

su 



(42) 



which means that normalization is conserved. For the 
second moment we obtain 



(x 2 (s,u)) = ^M[l + fi( u )] 



(43) 



which is the same as Eq. (|l7|). The fourth moment is 
more interesting: 

(x 4 ( Sl u)) = [1 + Sl(u)] \m 4 + 6fi(u)m?l . (44) 

u 

Higher order moments are obtained in a similar way, for 
the sake of space they are not included here. Odd mo- 
ments vanish due to the symmetry of the random walk. 

One can easily see that the ACTRW n th moment 
(x n (s,u)} depends on the microscopic jump moments 
{ni2, • • ■ , m n }. However in the limit u — > 



(x n (s,u)) 



h s (u) 



1 



u (Au a ) n / [ 



-.™r^ (45) 



which depends on mi but not on the higher order jump 
moments m^rriQ etc. Thus the moments m„ with n > 2 



where 

g(z) = z 
Hence 

P(x,t a ,t)\ x=0 



_ r q/2 Sm (gg) 



dyfe^. (48) 



27rr(l-a/2)y (1 + y)?/ 



2r(a)r(2-3a/2) \ t 



2r(l-a/2) 



(49) 

In the limit t >> t a we recover standard CTRW behavior 

In Fig. o, we present the behavior of the AC- 
TRW on the origin. The ratio of the scaled non- 
singular ACTRW Green function to the Montroll- Weiss 
non-equilibrium CTRW Green function on the origin, 
P(x,t,t a )t a / 2 \ x=0 /2r(l - a/2), is plotted versus the 
scaled dimensionless time t/t a . This ratio approaches 
one in the limit t >> t a , showing that the ACTRW 
process will converge to the standard non-equilibrium 
CTRW behavior when t >> t a . It can be clearly seen 
in the Fig. | that for a > 1/2, the ACTRW process 
has roughly converged to the non-equilibrium CTRW 
limit when t/t a ~ 1, while for a — * 0, the crossover 
to CTRW limit becomes extremely slow. For example, 
when a — 1/12, large deviations from the CTRW limit 
are clearly observed even when t/t a — 10 s . Since the 
limit a — > is important for several systems [|l7L [l5| , 
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it becomes clear that when a is small, the convergence 
towards the standard CTRW results becomes extremely 
slow, and the aging effect is of importance even when 
t > t a . 




FIG. 6: The behavior of the non-singular Green function 
on the origin normalized by the long time solution of the 
Montroll-Weiss non-equilibrium CTRW. The convergence of 
ACTRW towards the CTRW result is extremely slow when 
a -» 0. 

From Eq. (fl9|), we see that for f<(„ the non-singular 
part of P(x, t a , t)\ x= o increases with time t when a < 2/3. 
This unusual behavior is not unphysical, because the sin- 
gular delta function term is a decreasing function of time, 
and the total probability of finding the random walker in 
a small vicinity of the origin is decreasing monotonically 
with time as expected. 



V. BIASED ACTRW 

We now consider one dimensional biased ACTRW. We 
therefore use the small k expansion of f(k): 



f(k) = 1 + ikm\ ^" TO2 ' 



(50) 



where mi > is the averaged jump length. Differentiat- 
ing Eq. ([[2]) once with respect to k and taking k = we 
find the mean displacement of the random walker in s, u 
space: 



(x(s,u)) 



m\h B (u) 
u [1 — tp (u)] 



(51) 



where h s (u) is defined in Eq. ||). Differentiating Eq. 
( |l2| ) twice with respect to k, we find the second moment 
of the biased random walk 



(x 2 (s,u)} = 



h s (u) 



u [1 — ip (it)] 



2mi 



1 - V (u) 



m 2 



(52) 



1. Aging Einstein Relation 

We now derive a relation between the mean square dis- 
placement in the absence of bias, and the mean displace- 
ment of the particle in the presence of bias, reflecting 
the fluctuation-dissipation relation valid within linear re- 
sponse theory (see [l?], fl8| for related work) . The 
case of zero aging t a = were discussed in [Q, where 
some conceptual problems of linear response theory for 
systems exhibiting anomalous type of diffusion was dis- 
cussed (e.g., the non-stationarity of the process, the de- 
pendence of a on external field) . 

We assume that the random walk is on a one dimen- 
sional lattice with lattice spacing a, therefore 



fix) = P L S(x - a) + P R S(x + a). 



(53) 



1, hence the jump moments in Eq. 
-,2 



Here P L + P R 

are mi = (Pr — Pl)h and m-i = a 1 . We assume that the 
process obeys local detailed balance, namely Pl/Pr = 
exp (—aF/kbT) where T is the temperature. Using these 
conditions, and the assumption of weak field aF/kbT <C 
1, we have mi 
we find 



2 F/2k b T. Using Eqs. @ and @), 



{x [s 



F 
2k h T { 



(s,u))c 



(54) 



The subscript F in Eq. 



indicates the presence of 



external field F. (x 2 (s, u)) is the mean square displace- 
ment in the absence of a field, i.e. Eq. (E3) with mi = 
and m.2 = a 2 . Since the equation holds for the s,u do- 
main, it holds also for the t a ,t domain 



(x (£„)*)) J 



F 
2k h T 



[x 2 (t o ,t)) . 



(55) 



Thus the mean square displacement of the particle in 
the absence of the field (the fluctuation) yields the mean 
displacement in the presence of a weak field. When the 
waiting times are exponentially distributed, we obtain 
the usual Einstein relation between mobility and diffusion 
constant, which is independent of the age of the process 
t a ■ For experimental verification of Eq. ( 54 ) in the non- 
aging regime t a = and with a < 1, see [51, 52, |53"|. 



2. Asymptotic Behavior of Biased ACTRW 

From Eq. ( |5l| ) we can derive the behavior of the mean 
displacement in exactly the same way as done in Sec. 



EVA], and find 

(x (*«,*)> 



Ar(l+a) 



* > t a 

t < t a . 



(56) 



AT(a) 

For the second moment we use the small s, u behavior of 
Eq. p2 and find 



(x 2 (s, u)) 



1 



Au 1+a s a (u - s) 



2ml 
u a A 



m 2 



(57) 
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Inverting this equation using Eq. (|56j), we investigate 
now the dispersion: 



u (t a , t) ~ yj(x 2 (t a ,t)}-(x(t a ,t)} 2 . 



(58) 



Considering first the u <C s limit corresponding to t ^> 
f n we recover Shlesinger's result [Q 

cr (t„,t) ~ 



t?t 2a 
~A 2 ~ 



1 



_T(l + 2a) T 2 (l + a) 
hence if mi ^ and a < 1, one finds 
m 1 t a I 2 



m 2 i Q 



ATil + a)' 



<y (fa, t) 



1 



.4 



r(l + 2a) T 2 (l + a) 



(59) 
(60) 



For t ^> t a the dispersion of the biased CTRW grows like 
the mean Eq. (|56| ) , a behavior very different than normal 
Gaussian diffusion. 

Considering the s <C u limit of Eq. (|57| ) and using Eq. 
(pq), we find for aging limit t <C t a 



u (t a ,t) 



lm 2 tta' 1 2m 2 t^H 



l\t 2 t a 



AT (a) A 2 r(a)T(2 + a) A 2 T 2 (a) 

(61) 

Hence if mi ^ and a < 1, one finds 
o (t a ,t) ~ 



mi / 2 
^4~Y r(a)T(2 + a)2 



(62) 



which is valid for t, t a 3> A x l a and t < t a . As expected 
the dispersion decreases as age of the processes becomes 
older. Note that as a. — » 0, the first term in eq. ( |6l| ) 
becomes important. 
For a = 1 one finds 



(63) 



where in this case ^4 has the meaning of the mean time 
between jumps. The dispersion in this case is indepen- 
dent of the age of the system t a , as expected from normal 
diffusion. We see that the dispersion in normal diffusion 
processes is controlled by the second moment of jump 
lengths m 2 (even when mi ^ 0), while for CTRW and 
ACTRW mi ^ is the relevant parameter. 



VI. POSSIBLE APPLICATION: 
SCHER-MONTROLL TRANSPORT 

In this section we briefly point out to one possible ap- 
plication of ACTRW. Scher and MontroU B modeled 



transport in disordered medium based on CTRW the- 
ory. The fundamental reasons of why and when their 
approach is valid, while being the subject of much the- 
oretical research |35j ^6), are not totally solved. What 
is clear is that on a phenomcnological level, one can use 
the Scher-Montroll approach to fit behaviors of charge 
currents in a large number of experiments in very differ- 
ent systems. For example, transport of charge carriers 
in: organic photo-refractive glasses (57j], nano-crystalline 
Tj0 2 electrode s ]5q] , conjugated polymer system poly p- 
phenylene |5^, |30|1, and liquid crystalline xinc octakis Q . 

Scher and MontroU model such transport processes us- 
ing non- equilibrium biased CTRW theory with an effec- 
tive waiting time distribution. In experiments, this cor- 
responds to charge transport which is started at time 
t = 0, for example by a short photo flash applied on the 
system. After the initial triggering of the process, the 
charge carriers are transported using an external bias. 
For such initial conditions, we know that it is useful to 
assume that the physical transport process is described 
by the non-equilibrium biased CTRW. 

In an aging experiment one would start the process, by 
an external impulse (e.g., a photo flash), then wait for 
an aging period t a , and only after that period add the 
external bias. In this case, biased aging CTRW might 
become a useful tool describing the aging transport. At 
this time it is still an open question if ACTRWs can 
be used to describe aging in real systems. Further it 
is not clear if aging in the above mentioned systems |57], 
|58| , |59} p0| pl| is measurable, and if so do these very 
different systems exhibit any common aging effects in 
their transport? 

In the non-aging experiments (57], |58|, ^0|, ^l) one 
measures the current of charge carriers, which accord- 
ing to the predictions of Scher and MontroU exhibit a 
universal behavior: I(t) cx t a ~ 1 for short times and 
I(t) oc i - " -1 for long times. The transition time ti, be- 
tween these two behaviors depends among other things 
on the length of the system. The short time behavior 
corresponds to I(t) oc 4r(x), with the non-equilibrium 
CTRW behavior (x) oc t a , which yields immediately 
I(t) oc . The long time behavior is more compli- 

cated, and is due to absorbing boundary condition. 

In this paper we have calculated the mean displace- 
ment of the biased ACTRW, without including the influ- 
ence of the boundary. Thus we provide the aging cor- 
rections to the short time behavior of Scher-Montroll 
transport. According to Eq. (p36|), the Scher-Montroll 
behavior I(t) oc t a ~ 1 is replaced with I(t) oc t"^ 1 when 
f <C t a <C fr,- This behavior is independent of t, simi- 
lar to behavior of normal currents. In this aging regime, 
the current decreases as t a is increased, while in the non- 
aging case the current is decreasing when the forward 
time t is increased. The assumption made is that the 
dispersion of the probability packet during the aging pe- 
riod, is small compared with the length of the system. A 
detailed investigation of aging in Scher-Montroll trans- 
port systems will be given elsewhere. 
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VII. SUMMARY 

We have derived an exact expression for the Green 
function of ACTRW in Fourier-double Laplace space. 
This generalized Montroll - Weiss equation describes dy- 
namics of a large class of random walks. Since the CTRW 
describes a large class of physical and chemical systems, 
mainly disordered systems, we expect that ACTRW will 
be a valuable tool when aging effects are investigated 
in these systems. Interesting aging behaviors are found 
when the system turns non-ergodic, namely when the 
mean waiting time diverges, a < 1. We note that also 
when the mean waiting time is finite, aging behaviors 
may be observed, however only within a certain time win- 
dow. 

We showed that asymptotic behavior of the Green 
function is related to a few parameters of the underly- 
ing walk a, A, and m,2, while other informations con- 
tained in tp(t) and /(r) are irrelevant. The Green func- 
tion behavior is non Gaussian when a < 1, it is re- 
lated to Levy's generalized central limit theorem and 
to Gordeche-Luck's |H| fractal renewal theory. Unlike 
standard random walks or non-equilibrium CTRWs, the 
asymptotic Green function is a sum of two terms: a sin- 
gular term corresponding to random walks where number 
of jumps is zero and a non-singular term corresponding to 
random walks where number of jumps is one or more. Fi- 
nally, we note that the fractional Fokker-Planck equation 
framework ||, |6^], developed based on CTRW concepts 
p3| , can be modified based on the results obtained in this 
manuscript to include aging effects. This topic is left for 
a future publication. 



the random variable t^+i- It is given by: 
Pt a (tN+i) 



N+l 



= (<M * N+1 - J2 T A 1 & N < ta < *w+i))jv. (64) 



JV=0 



Here I(tjj < t a < tN+i) = 1 if tjy < t a < tjf+i, otherwise 
it is zero. The average in Eq. ( |64| ) is 



{■■■)n = (tiI + 1 1 I ^{n)dn---) 



(65) 



We consider the double Laplace transform t a — > s and 
tjv+i -> u of P ta (ijv+i) Eq. 0, P s (u). Using t N = 



Sili T i-> we nave 

e~ taS I (i N <t a < i N+1 ) dt a = 



g — tNS g — ijv + is g s 



(66) 



and 



-Utjv + l 



(- \ 

hl (5 t N+1 - t 4 dt N+ i = e u 2^i=i *. 

Using Eqs. @^), we have from Eq. @ 

OO 

Ps(u) = 



(67) 



JV=0 



Acknowledgments EB thank J. P. Bouchaud for point- 
ing out Ref. J3lj|. Discussions with L. Levitov and A. 
Hcucr motivated this work. 



VIII. APPENDIX A 

In this Appendix, we derive Eq. (||) using a method 
which is slightly different then the one used in |3l| . Con- 
sider a non-equilibrium renewal process which starts at 
time t = (in the ACTRW this initial time is —t a ). Let 
ii i = 1, ...N be dots on the time axis on which jumping 
events occur. Let ijv+i denote the time on which the 
first jump event occurred which is larger than t a , namely 
tjv < t a < Fjv+i. Note that N itself is not fixed, it is 
a random variable. The time intervals between jumps 
events are denoted by Tj = tj+i — U. 

The random variable we are interested in is t\, where 
ti = — t a . Since t a is a parameter in the problem, 
knowledge of statistical properties of Fjy+i yields the dis- 
tribution of t\. Hence let Pt a (fjv+i) denote the PDF of 



N+l 



(exp l-u n 



exp (-s £\ =1 nj - exp [~s X,i=i ' 



i=l 



(68) 

Using the fact that the random variables are indepen- 
dent and identically distributed, we have 

1 °° 

P s (u) = ~Y,bP N (u + s)tJj(u) - i, N+1 (u + a)] , (69) 



N=0 



where ib(u + s) = J °° exp[— (u + s)r]^i(r)dr. Summing 
Eq. © we find 



P s (u) 



ip(u) — ip(u + s) 1 



1 — ip(u + s) 



(70) 



Now the PDF ht a (ti) is obtained from Pt a (Fjv+i) us- 
ing t\ = fjv+i — < Q . According to definition of Laplace 
transform, we can write 

M*0 = 



11 



ir dtae ~ sta Jo 



dt a e- st « I dtie- u ^ht a (ti)>, (71) 

where C^ 1 is the double inverse Laplace transform. We 
use t\ = ijv+i — t a and find: 

M*i) = 



Hence it is easy to see that 



h s (u) = P s _„(u). 



(73) 



zr 1 



dt a e~ 



dhe- u ^-^P ta (iN+i) 
(72) 



Inserting Eq. © in Eq. @, we find Eq. (|). 
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